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Abstract 

We calculate numerically the electrical conductivity a, heat conductivity k and 
shear viscosity rj of the hot plasma present in the early universe for the temper- 
ature interval 1 MeV < T < 10 GeV. We use the Boltzmann collision equation 
to compute all the scattering matrix elements and regulate them by the thermal 
masses of the t- and u-channel particles. No leading order approximation is 
needed because of the numerical integration routines used. 
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1 Introduction 



The transport properties of the hot plasma present in the early universe are of great 
interest. The transport coefficients play a significant role in the phase transitions of the 
early universe [TJ, 0, |3|, M, the creation and development of the primordial magnetic 
fields ||, H and finally in the creation of the primordial density perturbations and 
therefore in the galaxy formation |7], ||. 

There has been many attempts to estimate the transport coefficients. (Textbook 
estimates for thermal conductivity, shear and bulk viscosity in the ultrarelativistic 
plasma were given in |J and |L0|.) In JO]] the viscosities of a pure gluon plasma and 
of a quark-gluon plasma were computed in the weak coupling limit from a variational 
solution to the Boltzmann equation. In JTIJ the transport coefficients were calculated 
for plasmas interacting through strong, electromagnetic and weak interactions to lead- 
ing order in the interaction strength. The rates of momentum and thermal relaxation, 
electrical conductivity and viscosities of quark-gluon and electro dynamic plasmas were 
also included in |12|]. In the transport coefficients and relaxation times were cal- 
culated to leading orders in the coupling constant for degenerate quark matter within 
perturbative QCD for temperatures and inverse screening lengths much smaller than 
the quark chemical potential. The viscosities of the QCD-plasma were considered in 
[Hjj . Textbook estimates for thermal conductivity and viscosity in the ultrarelativistic 
plasma were given in || and [K|. The conductivity of a relativistic plasma has also 
been considered in |H| by reformulating the collision operator for a relativistic plasma 
in terms of an expansion in spherical harmonics. 

However, all the previous estimates involve approximations. The calculations have 
been performed only to the leading (logarithmic) orders and certain scattering reac- 
tions present in the hot primordial plasma have been neglected. In the present paper, 
the only approximation we make is to assume that when particles appear in the heat 
bath of the primordial plasma, their thermal velocities are high enough to allow us to 
treat them ultrarelativistic and therefore massless. We consider radiation dominated 
plasma with temperatures from about 1 MeV to about 10 GeV. 

The viscous damping and heat conducting effects affect the first order phase tran- 
sitions in the early universe. During such phase transitions, instabilities may occur 
when the transport of latent heat is dominated by the fluid flow. In |J this was studied 
in the EW transition in the small velocity limit, in |fj] in the QCD transition, in || 
for cosmological detonation fronts and in [|J for general first order transitions with 
either large or small bubble wall velocities. The instabilities can be damped by finite 
viscosity and heat conductivity due to the diffusion of radiation on small length scales. 

When considering the creation and development of primordial magnetic fields it is 
also of importance to know the transport coefficients of the surrounding primordial 



plasma |5|, [J . Finite electrical conductivity leads to the diffusion of magnetic fields and 
therefore is of importance when trying to explain the further evolution of an initial seed 
magnetic field (see e.g. [16]). The instabilities in the first order phase transitions can 
be shown to create seed magnetic fields 0. Therefore the seed primordial magnetic 
fields also depend on the size of the transport coefficients in the primordial plasma in 
a crucial way. In particular, the seed fields created in the QCD phase transition are 
highly dependent on the neutrino viscosity ||. 

Galaxy formation |7| is likewise affected by the thermal coefficients in the primor- 
dial plasma || . Viscosity tends to heat up the plasma, and on the other hand thermal 
conduction transfers heat from regions of high temperature to regions at lower tem- 
perature. The effects affect the structure formation and thus in order to make precise 
models of the galaxy formation it is important to know the values of the transport 
coefficients. 

It should be noted that two major effects are involved when heat transport and 
viscosity are studied. These effects are the number of possible scattering reactions, 
i.e. the number of particles present in the plasma, and the Debye mass present in 
the propagator. As we will discuss in Section 4, the first effect will decrease the 
transport coefficients while the second effect will increase them. The net behaviour of 
the coefficients depend on the exact interplay between these two effects. 

The different interactions give contributions of different strengths to the electrical 



conductivity as was shown in |17| . However, in |L7j only strong and electromagnetic 
interactions were considered. Here we also deal with the scattering reactions involving 
neutrinos. 

The paper is organized as follows. In section 2 we introduce the reader to the 
method of calculation and, using the full definition of the energy-stress tensor T^ u , we 
define the transport coefficients presented here. In section 3 we reconsider electrical 
conductivity. We first discuss what is the proper way to introduce the Debye screen- 
ing to the propagators in the matrix elements. We then include additional scattering 
reactions not present in previous work and present a new value for the electrical con- 
ductivity. We also compare our results with other recent estimates for the electrical 
conductivity in the hot plasma. In Section 4 and 5 we apply our technique to calcu- 
lating thermal conductivity and shear viscosity in the hot plasma, respectively, and 
compare our results with other recent estimates. Section 6 contains a summary of the 
results and a brief discussion on their relevance for the early universe. 



2 The general formalism 



2.1 The energy-stress tensor 

We start by defining the transport coefficients. We do this by using the energy-stress 
tensor T^ u \ 

= (p + p)UW - + Tg + TZ, (1) 

where p is the total energy density, p the pressure, i7 M the velocity four-vector {U^U^ = 
— 1) and g^ v the metric tensor. is the usual electromagnetic energy-stress tensor 
(see e.g. |[L8|| ) and represents the non-ideal contributions from a finite viscosity 
and heat conductivity caused by the diffusion of the various particles present in the 
plasma. 

To find out the meaning of the transport coefficients we now follow || and derive 
the form of . Because of dissipation, we must also add a correction term to the 
particle current N^: 

= nU» + N£, (2) 

where n is the particle number density and N% the correction term. To avoid the am- 
biguities generated by the adding of extra terms induced by dissipation and diffusion, 
we should define what is meant by the number and energy densities: 

n = -EW„, 
p ee UJJpT* , 

ee {-N v N v y 1/2 N tl . (3) 

Inserting the definitions Eq. (^|) into Eq. (P and Eq. (0) we can easily see that the 
following conditions must be satisfied: 

K = o, 

U a U p Tf P = 0. (4) 
Using the conservation laws for the fluid: 

d u T^ u = , 



d v N v = (5) 



and making use of Eq. (pj) and Eq. (HI) we find that 



d v U" = -—d v n . (6) 



n 



Using Eq. (P, Eq. (||), Eq. (||) and the second law of thermodynamics, we can write a 
formula for the rate of entropy production using the total entropy current four-vector 
S^ 1 ee nsU^ 1 — U U T^/T, where s is entropy per particle and T the temperature: 

8^ = ~{d v U»)T£ + L {dvT )U,T^. (7) 



Now the task is to construct in the way that the rate of entropy production 
per unit volume given by Eq. (0) is positive for all possible fluid configurations. As 
explained in , we may take the dissipative term to be linearly dependent on the space- 
time derivatives of the four-velocity, energy and number densities. Also, the effects 
considered here are of first order and therefore the adiabatic equations of motion can 
be used to find out the form of . 

It is easiest to continue in the locally moving Lorentz- frame, where = (1, 0, 0, 0). 
As explained, can now be constructed as a linear combination of the derivatives 
dU^/dx", dU^/dt, dT/dx u and dT/dt. Of these, dU° jdx v vanishes in this frame and 
dT/dt can be expressed in terms of V-U because of adiabacity. Then the most general 
structure allowed by rotational and space-inversion is 

T tp = -vi^ + ^-ft-VsA -CV-U5*' , 



1 Tp 



(dU l 


dXJi 
dxi 


[dxj 


dT 


e- 


dxi 



T T ° P ° = 0. (8) 

By comparing Eq. with the usual non-relativistic hydrodynamics we can now recog- 
nize T), C and k to be the shear and bulk viscosity and the heat conduction coefficients, 
respectively. £ is a relativistic contribution with no non-relativistic counterpart. By 
using Eq. @ in Eq. fl7|) and demanding that the rate of entropy production is positive 
for all fluid configurations, we get the result £ = Tk together with the fact that all the 
transport coefficients are non-negative. 

Again, following ||, we consider a radiation dominated plasma having very short 
mean free times. The energy-stress tensor to first order in the free time r for such a 
fluid reads ||: 

' ^ u d P\ \ (TdW , rrY #T 



Tg= - \{^T*)/{^) n \[-— + W 



dT" ) \ 3 dx^i dx^ j 

/ 1 1 1 \ /^ r T 1 

AbTT 3 {[2U^U u U" 1 + -g» u W + -8^U U + -S^U" ' 



3 3 3 / <9x 7 

where b is a constant. Since is a tensor, it is sufficient to study it in the locally 
comoving Lorentz frame. Using once again the properties of this frame, U^ 1 = (1, 0, 0, 0) 
and dU° /dx^ = 0, we can somewhat simplify the energy-stress tensor, Eq. (|9]): 

n = -*TV[i-O„] 2 V.U-A6rvf^ + ^-^V.u)+0(r 2 ), 




TP 



= 0(r 2 ). 



(10) 



From Eq. (|T0D it is easy to extract the values for the transport coefficients: 

r) = (4/15)6T 4 r and 

k = (A/3)bT 3 r. (11) 

As can be seen from Eq. ([TT]), the transport coefficients are indeed proportional to 
the mean free time (or equivalently the mean free path) of particles in the plasma, in 
which interactions are sufficiently strong so that particles have a high scattering fre- 
quency. For various species of particles the coefficients are obviously different because 
their cross sections vary and thus they have different free paths in the plasma. In 
Sections 4 and 5 we calculate heat conductivity and shear viscosity considering all the 
particles and their scattering reactions in the hot primordial plasma. We note that 
the values obtained in Sections 4 and 5 are self-consistent with the assumption of the 
high scattering frequency, i.e. the mean free times estimated from the values of Table 
5 are very short. Let us finally note that because the plasma considered here is highly 
relativistic, and therefore p ~ |p, the bulk viscosity ( is approximately zero. 

2.2 The Boltzmann equation 

To actually calculate the transport coefficients, we need some tool to treat all the 
scattering reactions in the plasma. This tool is the Boltzmann equation, which in the 
expanding Robertson- Walker universe reads: 

f/ + p ' (I? - - e tH - § p2Rk - e w p0E ' = c(p - < )p0 • (12) 

where C (p, t) is the collision integral. Here we have assumed a Robert son- Walker 
metric with a scale factor R for the background, and we have defined Ff = —Ei 
and i™ = e y The Robertson- Walker metric is valid here if the mean free path 
of the particles is sufficiently small i.e. no large flows are present in the universe. 
The resulting values in Sections 4 and 5 suggest that the usage of this metric is a 
self-consistent assumption. We also prefer to use co-moving coordinates and define 

f(p,t)= J 5(p - (p 2 i? 2 + m 2 ) 1/2 )f(p,Po,t)dp . (13) 

Inserting this into Eq. (|T2|), making use of the local momentum defined as p = Rp 
and assuming adiabatic expansion in which E l ^> (R/R) |p| = H\p\ we can use the 



Lorentz-frame and write Eq. (|T^) in the form (from now on we drop the tildes and we 
also assume that the electric field is sufficiently small, E l <C T 2 ) 

To get a sufficiently good estimate, it is only necessary to consider 2— >2 reactions, and 
thus we may use as the collision integral the following expression: 

C(p,t) = — ( dP b dP c dP d (2ir) A 5 A (p + p b -p c -p d )\M(Ab -> cd)\ 2 F, (15) 
Po J 

where p is the four-momentum of the charged test particle A, and we have defined 
dPi = d 3 pi/ ((27i) 3 2E(pi)). The factor T contains the distribution functions of the 
initial and final states. The Pauli blocking factors must be included for the fermions 
and the enhancement factors for the bosons. Below are two examples of the jF-factor. 
The first is for the case when all the particles in the reaction are fermions, and in the 
second case particles b and c are bosons, such as is the case for Compton scattering: 

f [1 - /(p)][l - f b (p b )]fc(Pc)f d ( Pd ) - [1 - UPcW - fd(p d )]f(p)fb(p b ) , , , 

J Ui-/(p)][i + A(p*)]/^ 

The time dependence is not shown explicitly above. Let us now assume that there 
exists a constant electric field in the early universe (Of course, we do not claim such a 
field actually exists but rather use it as a probe of the plasma properties. Furthermore, 
we take it to be constant in the sense that its coherence length is larger than the mean 
free path of charged particles) and treat its effect as a perturbation on the distribution 
of the test particle. We write / = fo + Sf, where fo is the equilibrium distribution and 
Sf the small perturbation. Inserting this into the collision integral Eq. QT5D results in 

C(p,t) = -- f dP b dP c dP d (27r) 4 5 4 (p + p b -p c - Pd ) 
Po J 

x \M\ 2 ([f 0c f od [l-f ob ] + fo b [l-foc][l-f 0d ])Sf =C(p)Sf(p,t). (17) 

Eq. (|TT| ) assumes that all the particles in the reaction are fermions; generalization to 
other cases is straightforward. 

We have assumed that all the particles b, c and d have an equilibrium distribu- 
tion. Therefore, if also the test particle has an equilibrium distribution, the collision 
integral vanishes and only the perturbation term in Eq. ( |17D survives. Note that the 
equilibrium distribution depends only on the momentum. 

The Boltzmann equation can now be linearized in order to calculate the electrical 
conductivity. Assuming that Sf depends only on p, the linearized Boltzmann equation 
reads 

^- e W =C{p)sm < (18) 

where v 1 is the velocity of the probe particle a. Using Eq. (|1^) we can now begin to 
calculate the transport coefficients in the early universe. 



3 Eletrical conductivity 



There has been many attempts to estimate a. In (19j this was done in the relaxation 
time approximation and in |20[ and [^1 a Coulomb correction was applied. Recently a 
numerical work was carried out in |T7J which showed that the electrical conductivity in 
the early universe is mainly due to the leptonic contribution and its value for T < 100 
MeV a ~ 0.76T while at T ~ M w it was found that a ~ 6.7T. 



However, in |17fl the contribution from lepton-quark scatterings was neglected, 



which was correctly pointed out in p2| . We have now included these scatterings 
and thus acquired a better estimate for the electrical conductivity. The additional 
reactions now included are of the form I q — > l^q and I q — *> / ± g, and these affect only 
the leptonic contribution because the quarks interact via the strong interaction and 
thus the electromagnetic scattering processes mentioned above can still be neglected 
when considering the electrical conductivity created by the quarks. 

The t- and ^-channels give singular contributions to the collision integral and that 
is why one needs regularization. The full finite temperature calculations with higher 
order Feynman graphs and resummation would provide the regularization correctly 
but the expressions for the matrix elements become quite long. Another technical 
complication is the correct handling of the real and imaginary parts. A suitably 
accurate calculation can be, however, performed by using thermal masses as regulators 
in the t- and u-channel propagators. This approach can be expected to give roughly 
the same results as the one using the full thermal propagator. 

The thermal masses used in the t- and u-channels for the lepton, quark, photon 



and gluon propagators, respectively, are the following [23 



m ] = —T 2 ~ 0.0115T 2 , 



m 2 = 9s_ T i _ Q.251T 
q 6 

e 2 



m 2 = -(^Qf + 3Z q Q 2 )T 2 -OMOQ&Qf + 3£ 9 Q 2 )T 



3 

g 2 (3 + ^)T 2 ~ 1.508(3 + ^)T 2 . (19) 



where Nf is the number of quark families present, the sums are over all particles with 
m < T, and we adopt the values g 2 = 4na s ~ 1.508 and e 2 = 4na e ~ 0.0917, and 
T QCD = 200 MeV. 

The matrix elements for all the reactions used in the collision integral can be found 
in Tables 1, 2 and 3. All the elements are summed over final spins and averaged over 
initial spins; the QCD-processes are also summed over final color and averaged over 
initial color. Because we can treat every spin and color state as a different particle in 
the initial state we must, however, remove the averaging in the collision integral by 
multiplying the matrix elements with the number of initial spin and color states. 



The proper introduction of the thermal regulators in the propagators of the t- and 
w-channels should also be given special attention. For example, the matrix element for 
the reaction e~e + — *> e~e + is the sum (u 2 + s 2 )/t 2 + (u 2 + t 2 )/ s 2 + 2u 2 / (st), which in 
17] was first simplified into the form (s 4 + t A + u 2 (s + t) 2 )/(s 2 t 2 ) and only after that 



the regulators were inserted into the t- and w-channels. But one should be extremely 
careful when applying regularization because the adding of the regulators in different 
stages of the calculation changes the actual value of the matrix element, and therefore 
the value of the conductivity This is so because the regulators are added only to the 
t- and -u-propagators, not consistently everywhere in the matrix element. The proper 
way to introduce the screening terms is to apply the regulators to every term given by 
the Feynman graphs separately and before joining the terms together. The regulating 
should be done this way because the different terms in the matrix elements have 
singularities of different order and hence the terms should be regulated and integrated 
separately. 

We consider the temperature interval from about 1 MeV to about 10 GeV and make 
the simple approximation that particles appear in the thermal bath when temperature 
is greater than their mass. The collision integral is performed numerically by evaluating 
the integral by a simple Monte Carlo simulation. 

The results, however, are not very different from the previous ones ||17|| ; the quark 
contribution is still negligible compared with the leptonic electrical conductivity and 
the actual value of the conductivity is still ~ 10T in natural units. The effect of 
neutrino scatterings to the contibution of either QED or QCD electrical conductivity 
is negligible. The results are given in Fig. 1 and in Table 5. 

The main difference between the value for the electrical conductivity here and 
in p2[ is due to the difference in the definition of the electrical conductivity. In 



2^1 an electric field extending itself over the whole universe is used to probe the 



electrical conductivity of the early universe. Because in p^] a homogenous field is 
assumed to cover the entire space, the leptons (quarks) and antileptons (antiquarks) 
are accelerated in the opposite directions. In this scenario, it is justified to neglect 
in the first order all other scattering reactions than the lepton-antilepton (quark- 
antiquark) annihilation. But the resulting electrical conductivity is a global one - one 
that is the same in the whole universe and does not take into consideration the small 
patch-like structure of the primordial magnetic field. 

The primordial magnetic field B, if it exists, is most likely random: both its mag- 
nitude and direction vary. Therefore in the early universe it is more sensible to talk 
about the local electrical conductivity, 0"(x). The magneto hydrodynamic equation 
reads then 

f. V x ( vxB)-Vxfe?). m 



From Eq. p0| ) we see that to blindly treat the electrical conductivity strictly as a con- 
stant would disregard some aspects of the dissipation of the magnetic field. However, 
to a good approximation we can take the electrical conductivity to be a constant inside 
the patches of the magnetic field characterized by correlation length L. 

The magnetic field is coherent inside the correlation length L. The different patches 
containing magnetic field are uncorrelated and thus a convenient way to define the local 
conductivity is the measure of the dissipation of the magnetic fields in the coherent 
patches. We measure the local electrical conductivity as experienced by a charged 
particle shot out of one of the magnetic patches to the next, uncorrelated one. The 
particles are ultrarelativistic because of the high temperature and therefore the differ- 
ent bulk velocities between the magnetic patches can be neglected. 

The condition when this scenario is valid can be easily checked. The mean free path 
of a particle in the plasma is lf ree = (cr ip n p ) _1 , where a tp is the transport cross section 
of the particle species p, and n p is its number density. The condition for validity is 
now that the mean free path is longer than the coherence length of the magnetic field, 
L < lf r ee- Typical coherence length could be L ~ 1/T whereas lf ree — l/(a 2 T). In 
this case particles will typically be shot out of a given patch before they have time to 
interact. When they finally do, this will happen in a different patch where the other 
particles have a bulk motion of their own. 

However, the coherence length of the primordial magnetic field is not known exactly. 
It might very well be of the order of the natural microphysical scale, the interparticle 
distance, but there could be processes and phenomena that would cause it to be much 



larger. For example, one such process studied recently is the inverse cascade [p4 
which transfers magnetic energy from small length scales to larger length scales. 

Here we shall assume that the mean free path of a particle is longer than the 
correlation length of the primordial magnetic field. To calculate the local conductivity, 
we assume one test particle from each species of particles being shot out of one patch 
of magnetic field to another, totally uncorrelated one. In the new patch, we assume a 
patch-wide electric field to give us a tool to measure the conductivity. Of course, this 
is just a test field, we do not claim such fields to exist in the early universe. Inside the 
patch, we also assume isotropy and therefore we can write Eq. fllTf) to the form 



e w El = d{p)5f{p) - (21) 



From Eq. (BTf) one easily obtains 



*/(p) = -^^tf. (22) 
v ; C(p) dp 1 v ' 



The induced current density is given by 

5fip)P d* P , (23) 
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Figure 1: a/T as a function of temperature, 
and conductivity a a, associated with a given particle species A, is defined by 

U = ^aE . (24) 



Thus the contribution of a single species A to conductivity is a A ~ 1/ \M(AX — > 
Y)\ 2 , where the sum is over all the processes which scatter the test particle A. For the 
purpose of conductivity, we may view the mixture of different particle species of the 
early universe a multicomponent fluid. The flow of each component contributes to the 
total current and adds up to the total conductivity, which reads 

vtot = a A , (25) 

A 

where the sum is over all the relativistic charged species present in the thermal bath. 
Note that the total conductivity is dominated by the species that has the weakest 
interaction. This is because the weaker the interaction, the longer time the current 
flow is maintained. 

We consider the temperature interval 1 MeV < T < 10 GeV, and make the simple, 
crude assumption that particles appear in the thermal bath only when temperature 
is greater than their mass. Thus below 100 MeV, for example, the only charged 
particles present are the electrons and positrons. When T > Tqcd, also the quarks 
should be counted in. Their main interactions are strong, so that their electromagnetic 
interactions may be neglected. The list of the relevant reactions involving QED and 
QCD charged particles can be found in Tables 1 and 2. 



4 Thermal conductivity 



To find out the thermal conductivity, we now abandon the isotropy requirement and 
the requirement of locality used in the previous Section in computing the electrical 
conductivity Let us for a while let the equilibrium distribution f depend also on the 
chemical potential //: /o = 1/(1 ± e^ _M ^ T ). The conditions for thermal equilibrium 
require the constancy of temperature and of the sum fi+V, where V is the energy of the 
particles in an external field, throughout the media considered. Here we take V = e<f>, 
where <fi is the electric field potential. In a plasma with a non-uniform temperature 
distribution, the electric field E is not zero even if the current is zero. In general, when 
both the current density j and the temperature gradient VT are not zero, the relation 
between these quantities and the electric field can be written as 

j = aE-aaVT+ (26) 

e 

where a is the electrical conductivity and a the thermoelectric coefficient. We can 
now write the Boltzmann transport equation in the form 

-eE.^ + V.^ = C(p,T)5f. (27) 
ap p or 

Assuming that the temperature T and the chemical potential fi are functions of the 
spatial coordinates it is easy to solve 5f, the small perturbation to the equilibrium 
distribution / from Eq. ( pTf) : 

Sf = -{-eE^ + f^ + ^VT} 

pC(p,T) op dfi dT 

± e (p-n)/T sjt 
= 7 7 ^ {eE + Vu+(p-u) } (28) 

In order to determine the thermoelectric coefficient a, we must now assume for a 
while that the condition eE + V/i = holds. From the definition of the electric 
current Eq. (^3|) and from Eq. (|26|) one gets 

j = -aaVT = -e [ ^d 3 p. (29) 
J p 

Remembering the earlier definition of a, Eq 
the thermoelectric coefficient: 

a = 

where 



G = 



1 ± e(p-ri/ T ) 2 C(p) 





(32) 



From the energy flux relation 



q = / p5fd 3 p EE -KVT, 



(33) 



and from noting that in the early universe we can put the chemical potential jj to zero, 
we get an expression for the thermal conductivity coefficient k: 



The plus-minus sign refers to fermions and bosons, respectively. The heat conductivity 



its sign in the collision integral when considering fermions or bosons, respectively 

As explained in Section 2, thermal conductivity is related to the mean free time, or 
equivalently, to the mean free path. The values for the heat conductivity obtained in 
this Section are self-consistent with the assumption of short mean free times in Section 
2 as can be seen by comparing the values presented in Table 5 and Eq. (|TT|) . 

From Eq. (|34] ) it can be easily seen that k is proportional to the inverse of the square 
of the coupling constant of the scattering reaction in question. Thus the smaller the 
coupling constant is the bigger the value for k will be i.e. k indeed is related to the 
mean free path. 

The subtle interplay between the increasing effect on thermal conductivity by the 
temperature dependent regulators and the decreasing effect due to the change in the 
amount of particles (potential scatterers) present should also be noted. As can be 
seen in Fig. 2 and Table 5, the case for the leptonic thermal conductivity, Ki, is clear. 
The first effect dominates in the lepton case, increasing ni/T 2 all the way from the 
case where only electrons are present to the case where also the b-quark is present 
in the heat bath. The photon case is also straightforward. Because there are no 
temperature dependent regulators the only effect comes from the fact that more and 
more scatterers are present the higher the temperature is. Thus we see in Fig. 2 and 
Table 5 that the photonic thermal conductivity divided by the temperature squared, 
k 7 /T 2 , is decreasing steadily as the temperature rises. 

The case is similar with gluons and quarks as can be seen in Fig. 3 and Table 
5. The different behaviour of the quark thermal conductivity K q and the gluonic 
thermal conductivity n g can be understood when one looks at the different processes 
that are involved in the scatterings. All QCD matrix elements are listed in Table 2. 
The quarks get the biggest contribution to their collision integral from the potentially 




(34) 



coefficient is, however, always positive because the T (Eq. flTE|)) factor also changes 
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Figure 2: k/T 2 as a function of temperature for leptons and photons. 

singular matrix element which has been regulated by thermal masses while the gluons 
are not that sensitive to thermal masses. 

The thermoelectric coefficient a in Eq. ( PH| ) for leptons is found to be ±39.1 (de- 
pending on the sign of the lepton electric charge) and for quarks 47.4/e g , where e q is 
the electric charge of the quark. 

The neutrinos do not have any singularities in their matrix elements and thus as 
the temperature grows and more and more scatterers appear in the heat bath the 
neutrino thermal conductivity divided by the temperature squared, k u /T 2 , decreases 
steadily k u /T 2 is shown in Fig. 4. 

Earlier attempts to estimate the heat conductivity can be found, for example, in 
[ini II2L Ol. However, in fllOl no actual matrix elements were calculated but heat 



conductivity was estimated to be proportional to the mean free time of a particle 



species studied. In [O] and a more thorough calculation was performed to the 
leading logarithmic order. However, due to the approximations used in the formalisms 



of [ 1 2 1 and [ 1 3 1 we believe that values in the present paper are more accurate. 



5 Viscosity 

To find out the shear viscosity in the early universe we, once again, set the stage by 
starting from the Boltzmann equation Eq. (^) and the definition of the shear viscosity 
and assuming stationary viscous flow: 



v(diV k + d k Vi - j^ fc V • V) = J d 3 PPt v k 6f 



(35) 
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Figure 3: k/T 2 as a function of temperature for quarks and gluons. 
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Figure 4: k/T 2 as a function of temperature for neutrinos. 



Table 1: Matrix elements |M| 2 /(327r 2 a eTO ) for the QED processes used in the collision 
integral. All matrix elements are summed over final spins and averaged over initial 
spins. a em is the electromagnetic coupling constant, i ^ j, and e q is the charge of the 
quark in question. 
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Table 2: Matrix elements |M| 2 /(167r 2 a s ) for the QCD processes |25] used in the colli- 
sion integral. All matrix elements are summed over final spins and colors and averaged 
over initial spins and colors. a s is the strong coupling constant, and i ^ j. 
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Table 3: Matrix elements \M\ 2 /G 2 for the neutrino processes used in the collision 
integral. All matrix elements are summed over final spins and averaged over initial 
spins. G is the Fermi coupling constant, and i ^ j. For other abbreviations, see Table 
4. 
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Table 4: Abbreviations for Table 3 
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where i] is the shear viscosity, v the velocity of a particle and V the velocity of the 
flow. As in Section 4, we do not need the requirement of locality used in Section 3 in 
the computation of the electrical conductivity. 5f can now be deduced from Eq. ([Tj|) 
with E = 0: 

*/=^Ll (36) 
C(p) 

Substituting Eq. (KB) into Eq. (R3f) we see that 



/ rfW,^f ^ - / ^Wv • V(p • V)^-^- (37) 
J C(p) J dp C(p) 

Let us now assume a flow V = aye x , where a is a constant. Using V, Eq. (|35| ) and 
Eq. (|37|) we can write the following expression for the shear viscosity: 

d 3 pv x v v p x p v — — — — . (38) 

Using the fact that v and p have the same directions and therefore piVk = PkVi, and 
performing the integration in spherical coordinates, we get 



V 



] /*7T r1-K roo -|-pP/T 

- de d<$> p 4 sin 2 $ cos 2 $ sin 5 9-; . (39) 

TJo Jo Jo F C(p){l±eP/ T ) 2 



The ± sign stands for fermions and bosons, respectively. The shear viscosity is, how- 
ever, always found to be positive. Just as in the case of the thermal conductivity, 
we see here that the smaller the interaction of the particle species in question is, the 
bigger the shear viscosity of that species in the plasma is i.e. the viscosity is related to 
the mean free time of the particles, too. Here we see a similar behaviour as with the 
thermal conductivity in Section 4: the viscosity depends on the number of scattering 
reactions and on the thermal regulators. The scattering reactions are the same as in 
the previous sections, presented in Tables 1, 2 and 3. In Fig. 6 we show the results for 
the shear viscosity of the leptons and photons, and in Fig. 7 the shear viscosity of the 
quarks and gluons are presented. In Fig. 8 the neutrino shear viscosity is presented; 
all the results are collected in Table 5. It should be noted that the shear viscosities are 
the same for particles and antiparticles. Quarks and gluons interact strongly and thus 
their shear viscosity is the smallest while, respectively, the neutrinos interact weakly 
so that their shear viscosity is the largest. This can be understood by considering 
the shear viscosity with the following example. Let us shoot a particle into the hot 
(primordial) plasma. The weaker its interactions with the surrounding plasma are, the 
longer it will travel without scatterings, and vice versa. This is consistent with how 
we defined the shear viscosity in Section 2: it is indeed related to the mean free time 
(the mean free path) of a particle in the plasma. The values for the shear viscosity 
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Figure 5: rj/T 3 as a function of temperature for charged leptons and photons 

obtained in this Section are self-consistent with the assumption of short mean free 
times in Section 2 as can be seen by comparing the values presented in Table 5 and 
Eq. (0). 

Earlier estimates for the shear viscosity in hot (primordial) plasma can be found, 
for example, in || 11], [L2|, [13], Q . However, in || no matrix elements were calculated, 
and it was only noted that the shear viscosity is proportional to the mean free path of 
a particle species studied. O, O, O, H contained much more thorough calculations 



of the shear viscosity to the leading logarithmic order. The results, however, differ 
somewhat from the ones presented in Table 5 and Figs. 5, 6 and 7. We believe 
that the proper introduction and numerical integration of the matrix elements in the 
collision integral in the present calculation yields more accurate values for r\ than in 
the earlier publications. 

As noted before, in the early universe p ~ |p, so that the bulk viscosity (Eq. (|TT])) 
is approximately zero. 



6 Summary 

We find that the electrical conductivity is dominated by the leptonic current and that 
the quark and neutrino contribution to the electrical conductivity can be neglected. 
The values for the electrical conductivity can be found in Table 5. The definition of the 
electrical conductivity should be carefully considered. We started from the magneto 
hydrodynamic equation Eq. (|20|) and noted that the electrical conductivity defines the 
dissipation rate of the magnetic field in the plasma. But because of the short corre- 
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Figure 6: r]/T 3 as a function of temperature for quarks and gluons 
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Figure 7: r]/T 3 as a function of temperature for neutrinos. 



lation lengths of the magnetic fields in the early universe, the electrical conductivity 
must be a local quantity. The different patches of the magnetic field are uncorrelated 
and therefore our way of defining and calculating the electrical conductivity seems 
natural: when the mean free path of the particles in the hot primordial plasma is 
longer than the correlation length of the magnetic field, particles will be shot out of 
the small magnetic patches to the nearby uncorrelated ones and we can calculate the 
electrical conductivity by examining how fast a charged particle loses its initial bulk 
momentum. To define a global electric field which would extend itself over the whole 
universe would disregard the small random structure of the primordial magnetic field. 
We assume in the paper that the mean free path of the particles studied in the primor- 
dial plasma is always longer than the correlation length of the primordial magnetic 
field. 

To compute the heat conductivity and shear viscosity we use the same tools as in 
the computation of the electrical conductivity, but without the requirement of locality. 
The creation of the primordial magnetic fields can also depend on the viscosity and 
heat transportation of the hot primordial plasma || || . The viscous damping and heat 
conducting effects do also affect the first order phase transitions in the early universe 

[H] . During the phase transitions, instabilities may occur when the transport of latent 
heat is dominated by the fluid flow. The instabilities can be damped by finite viscosity 
and heat conductivity due to the diffusion of radiation on small length scales. The 
galaxy formation can also be considered to be dependent on the transport coefficients 
in the early universe [|7|, ||. A summary of the values for the heat conductivity and 
shear viscosity is given in Table 5. 

Reliable estimates for the transport coefficients will help us to understand better 
the creation of the primordial magnetic fields, phase transitions in the early universe 
and the creation of density perturbations in the primordial plasma, which then would 
seed galaxy formation. 
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